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Abstract 

We introduce an algorithm that embeds a given 3-connected planar graph as a convex 
3-polytope with integer coordinates. The size of the coordinates is bounded by 0(2^'^^") = 
0(188"). If the graph contains a triangle we can bound the integer coordinates by 0(2'* *^"). 
If the graph contains a quadrilateral we can bound the integer coordinates by 0{2^'^^^). The 
crucial part of the algorithm is to find a convex plane embedding whose edges can be weighted 
such that the sum of the weighted edges, seen as vectors, cancel at every point. It is well 
known that this can be guaranteed for the interior vertices by applying a technique of Tutte. 
We show how to extend Tutte's ideas to construct a plane embedding where the weighted 
vector sums cancel also on the vertices of the boundary face. 

1 Introduction 

Problem Setting. The graph of a polytope is an abstraction from its geometric realization. 
For a 3-polytope, the graph determines the complete combinatorial structure. The graphs of 3- 
polytopes are characterized by Steinitz' seminal theorem |31j . which asserts that they are exactly 
the planar 3-connected graphs. 

A natural question is to ask for a geometric realization of a 3-polytope when its combinatorial 
structure is given. One might be interested in a realization that fulfills additionally certain op- 
timality criteria. For example a good resolution is desirable to obtain aesthetic drawings [HJ 130) . 
We address a different problem and ask for an embedding whose vertices can be placed on a small 
integer grid. The vertex coordinates of such an embedding can be stored efficiently. 

Related Work. Suppose we are given the combinatorial structure of a 3-polytope by a graph G 
with n vertices. The original proof of Steinitz' theorem transforms the 3-connected planar graph G 
into the graph of the tetrahedron by a sequence of elementary operations. The transformation pre- 
serves the realizability as a 3-polytope. Since all operations can be carried out in the rationals, the 
proof gives a method to construct a realization of a 3-polytope with integer coordinates. However, 
it is not easy to keep track of the size and the denominators of the coordinates, which makes it dif- 
ficult to apply this approach for our problem. An alternative proof of Steinitz' theorem goes back 
to the Koebe-Andreev-Thurston Circle Packing Theorem (see for example Schramm f55^). This 
approach relies on non-linear methods, which make the (grid) size of the embedding intractable. A 
third proof of Steinitz' theorem relies on the "liftability" of planar barycentric embeddings. Since 
this barycentric approach is based on linear methods, its construction favors computational aspects 
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of the embedding. This led to a series of embedding algorithms: Hopcroft and Kahn [T^, Onn and 
Sturmfels [20^, Eades and Garvan |12J, Richter-Gebert (24j, Chrobak, Goodrich, and Tamassia [3]. 
Our work also follows this paradigm. 

As a first quantitative analysis of Steinitz' theorem, Onn and Sturmfels |20) showed that integer 
coordinates smaller than n^^^" suffice to realize a 3-polytope. Richter-Gebert improved this bound 
toO(2i8" ). A more careful analysis of Richter-Gebert's approach shows that the size of the integer 
coordinates can be bounded by 2^^" |21) . 

Integer realizations with at most exponentially large coordinates in terms of n were previously 
known for polytopes whose graph contains a triangle (Richter-Gebert [21]). We describe this 
method in Section 3.1 (p.jsj) as Case 1 of our embedding algorithm. In Richter-Gebert's approach 
(and already in Onn and Sturmfels [20]), graphs without triangles are embedded by first embedding 
the polar polytope, whose graph in this case has to contain a triangle. Based on the polar, an 
embedding of the original polytope is constructed. However, this operation yields coordinates with 
a quadratic term in the exponent. 

For triangulated 3-polytopes, Das and Goodrich |10j showed that they can be embedded with 
coordinates of size , using an incremental method which can be carried out in 0{n) arithmetic 
operations. Triangulated 3-polytopes are easier to realize on the grid than general polytopes, since 
each vertex can be perturbed within some small neighborhood while maintaining the combinatorial 
structure of the polytope. An explicit bound on the coordinates has not been worked out by the 
authors. For stacked polytopes a better upper bound exists [,38) . but it is still exponential. 



Lower Bounds. Little is known about the lower bound of a grid embedding of a 3-polytope. An 
integral convex embedding of an n-gon in the plane needs an area of fJ(n^) [U [S] |32l [35] . Therefore, 
realizing a 3-polytope with an (n — l)-gonal face requires at least one dimension of size r2(n'^/^). 



Two Dimensions. In the plane, planar 3-connected graphs can be embedded on a very small 
grid. For a crossing-free straight-line embedding an 0{n) x 0{n) grid is sufficient |l^[27j. This is 
also true if the embedding has to be convex A strictly convex drawing can be realized on an 
0{n^)x 0(n2) grid [3 . 



Higher Dimensions. Already in dimension 4, there are polytopes that cannot be realized with 
rational coordinates, and a 4-polytope that can be realized on the grid might require coordinates 
that are doubly exponential in the number of its vertices. Moreover, it is NP-hard to even decide 
if a lattice is a face lattice of a 4-polytope [211 [2S] . 



Results. In this article we develop an algorithm that realizes G as a 3-polytope with integer 
coordinates not greater than 0(187.13") — 0(2^-^'^"). This implies that for any 3-polytope a 
combinatorially equivalent polytope can be stored with 0{n) bits per vertex. For the case that G 
contains a triangle we show that G admits an integer realization with no coordinate larger than 
0(28.4") = 0(2*-*^"), if G contains a quadrilateral face, the size of the coordinates can be bounded 
by 0(43.99") = 0(2^-^'^"). The most difficult part of the algorithm is to locate the boundary face 
of the plane embedding such that a lifting into exists. This problem can be reduced to a 
non-linear system which is most complex when G contains neither a triangle nor a quadrilateral 
face. 

Partial results containing the essential ideas for graphs with quadrilateral faces (Case 2 of Sec- 
tion 3.1 1 were presented by the second author at the workshop The Future of Discrete Mathematics 



at Stifm Castle, Czech Republic, in May 1997. The results of this paper were presented in a dif- 
ferent form at the 23rd Annual Symposium on Computational Geometry in Gyeongju, Korea, in 
June 2007 |22j . Since then, we were able to simplify the computation of the explicit bounds with 
help of Lemma |3.10[ The simplification yields slightly different bounds. By improving the bound 
of Lemma 9 in |22j by a polynomial factor (now Lemma 3.9 ) we obtain better bounds in the end. 
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However, our analysis could be further improved with help of the more complicated construction 
of [22j. Since the improvement would only result in a constant factor we decided to present the 
simpler and more elegant analysis. 

A follow-up work [30\ extends the techniques of this article and studies more general barycentric 
embeddings. With help of these modifications, a grid embedding with z-coordinates smaller than 
2n can be constructed. The small x-coordinates are realized at the expense of the size of the y 
and z-coordinates, which are bounded by 2'-"-^ log"). 

Remark: Most recently, Buchin and Schulz [5] improved the upper bound for the maximum 
number of spanning trees contained in a planar graph. This has a direct consequence for our results, 
since we obtain the bound for the necessary grid size in terms of this quantity. In particular, the 
new bounds of [5] yield that our algorithm requires a grid of size 0(147.71") = 0(2^-^^") (general 



case), 0(39.87") = 0(25-32") (g- contains a quadrilateral face), and 0(27.94") = 0(24-8i") (G 



contains a triangular face). 

2 Lifting Planar Graphs 

Let G — {V, E) be a 3-connected planar graph with vertex set V — {wi, . . . , w„} embedded in the 
plane with straight edges and no crossings. The coordinates of a vertex Vi in the (plane) embedding 
are called Pi := {xi,yi)'^, the whole embedding is denoted as G(p). Let h: V ^ M.he a height 
assignment for the vertices in G. We write Zi for h(vi). If the vertices {xi,yi, Zj) of every face of 
G lie on a common plane, we call the height assignment h a lifting of G(p). 

Definition 1 (Equilibrium, Stress). An assignment w: i? — )■ M of scalars {denoted as uj(i,j) = 
ujij — ujji) to the edges of G is called a stress. 

1. A vertex Vi is in equilibrium in G(p), if 



2. The embedding G(p) is in equilibrium if all vertices are in equilibrium. 

3. If G(j3) is in equilibrium for the stress uj, then uj is called an equilibrium stress for G(p). 

It is well known that equilibrium stresses and liftings are related. Maxwell observed in the 19th 
century that there is a correspondence between embeddings with equilibrium stress and projections 
of 3-dimensional polytopes [T^. There are different versions of Maxwell's theorem. For the scope 
of this article the following formulation is the most suitable. 

Theorem 2.1 (Maxwell, Whiteley). Let G be a planar 3-connected graph with embedding G(p) 
and designated face fi . There exists a correspondence between 

A) equilibrium stresses uj on G(p), 

B) liftings o/G(p) in M.^ , where face fi lies in the xy-plane. 

The proof that A induces B (which is the important direction for our purpose) is due to Walter 
Whiteley (5B] . The Maxwell-Cremona correspondence finds interesting applications in different 
areas (see for example Hopcroft and Kahn [14], and Connelly, Demaine and Rote ^). 

To describe a lifting, we have to specify for each face fi of the graph the plane Hi on which it 
lies. We define Hi by the two parameters a, and di. The plane Hi is characterized by the function 
that assigns to every point p in the plane a third coordinate by 





(1) 



iJj! p 1-^- (p,a.,) +di. 



(2) 
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Here, (•, •) denotes the dot product. The correspondence between hftings and stresses comes 
from the observation that the "slope difference" a; — between two adjacent faces // and fr is 
perpendicular to the edge — pj that separates them: 

ai - Sir ^ ujij{pi - pj)^, (3) 

for some scalar ojij S E. Here, p""^ :— (^^) denotes the vector p = rotated by 90 degrees. It is 
not hard to show that these numbers form an equilibrium stress. 

The other direction, the computation of the lifting of G{p) induced by uj is straightforward, 
see Crapo and Whiteley We follow the presentation of Connelly, Demaine and Rote [7 for the 
computation of the lifting. 

The parameters a^ and di can be computed by the following iterative method: We pick fi as 
the face that lies in the xy-plane, and set ai = ([]) and do = 0. Then we lift the remaining faces 
one by one. This is achieved by selecting a face /; that is incident to an already lifted face J^. Let 
be the common edge of /; and fr- Assume that in G{p) the face /; lies left of the directed 
edge ij, and fr lies right of it. The parameters of Hi can be computed by 

ai ^ ujijipt - pj)^ + ar, (4) 
di = ujij{pi,pf) +dr. (5) 

The formula Q comes directly from (|3|, and ([5| comes from the fact that the two planes must 
intersect above p^ and Pj. 

The sign of the stresses allows us to say something about the curvature of the lifted graph. 
According to ([3| and Q , the sign of ojij that separates fi and fr tells us if the lifted face fi lies 
below or above Hr. As a consequence we obtain the following: 

Proposition 1. Let G{p) be a straight-line embedding of a planar 3-connected graph G with 
equilibrium stress. If the stresses on the boundary edges are negative and all other stresses positive 
then the lifting induced by such equilibrium stress results in a convex 3-polytope. 

Lemma 2.1. // G(p) has integer coordinates only and the equilibrium stress is integral on all 
interior edges, then the z-coordinates of the lifted embedding are also integers. 

Proof. We select an interior face as face /i. The gradient ai — (0,0)"^ and the scalar di = 
are clearly integral. For all other interior faces fi the parameters a^, di of the planes Hi can 
be computed with help of equations Q and ([S]). By an inductive argument these parameters 
are integral as well. Computing the z-coordinate of some point Pi by ^ boils down to the 
multiplication and addition of integers. □ 



3 The Grid Embedding Algorithm 
3.1 The Plane Embedding 

The embedding of G as a 3-polytope uses the following high level approach. First we embed G 
in the plane, such that it is liftable (see Section pj), then wc lift the embedding to M?, finally 



we scale to obtain integer coordinates as described in Section 3.2 The analysis of the algorithm 
in Section |3.3| gives the new upper bound. The most challenging part is to construct a liftable 
2d-embedding. 

An embedding is called barycentric if every vertex that is not on the outer face is in the barycen- 
ter of its neighbors. Tutte showed that for planar 3-connected graphs the barycentric embedding 
for a fixed convex outer face is unique |33l I34j . Moreover, if embedded with straight lines, no two 
edges cross, and all faces are realized as convex polygons. In the barycentric embedding all vertices 
that are not on the outer face are in equilibrium according to the stress w = 1. Our embedding 
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Figure 1: A small example graph. 



algorithm uses this special stress only, although we state the lemmas as general as possible. Since 
our techniques might find applications in other settings we develop our main tools for arbitrary 
stresses. (Note that Tutte's approach works with arbitrary stresses that are positive on interior 
edges, see for example Gortler et al. [13].) 

We describe now how to compute the barycentric embedding of G. Let /o be a face of G that 
we picked as the outer face, and let k be the number of vertices in /q. For simplicity we want k 
as small as possible. Euler's formula implies that every planar and 3-connected graph has a face 
/o with fc < 5 edges. We assume that the vertices in G are labeled such that the first k vertices 
belong to /q in cyclic order. Let B :— {1, . . . , fc} be the index set of the boundary vertices and 
let / := {fc + 1, . . . , n} denote the index set of the interior vertices. The edges of /o are called 
boundary edges, all other edges interior edges. The stresses on the exterior edges will be defined 
later, but since they don't matter for the barycentric embedding we set them to zero for now. 

We denote with L = (lij) the Laplacian matrix of G (short Laplacian), which is defined as 
follows 

' if (i, j) e E and i ^ j, 

otherwise. 



1 the Laplacian equals the negative adjacency matrix of G.with 
We subdivide L into block matrices indexed by the sets / and B, 
and obtain Ljg, Lgj, Lbbi and Ljj. The matrix Ljj is called the reduced Laplacian matrix of G. 
For convenience we write L instead of L77 



For the special "weights" u = 
vertex degrees on the diagonal. 



Example. Consider the graph of Figurejl] with B — {1, 2, 3} and / = {4, 5, 6, 7}. We have 



/ 2 
2 
2 



L = 



-10 
-1 -1 
0-1 
V 0-1-1 



-1 -1 0\ 
0-1 0-1 
0-1-1 

■3-1-1 

-1 4 0-1 
-1 3-1 
0-1-1 4/ 



Lbb Lbi 
Lib Ljj 



with L = Ljj = 



I 3-1-1 0\ 
-1 4 0-1 
-1 3-1 
V -1 -1 4 / 



In the example the presence of the boundary edges (1,2), (1,3), and (2,3) is not refiected in the 
Laplacian, because the stress is set to zero on the boundary. 

The location of the boundary vertices is given by the vectors = (xi, . . . , zj.)"^ and = 
(yi, . . . , j/fc)"^- Since every vertex should lie at the (weighted) barycenter of its neighbors, the 
coordinates of the interior vertices x/ = (x^+i, . . . , a;„)^ and yj = (?/fe+i, . . . , yn)^ have to satisfy 
the equilibrium condition ([T]) for the stress uj. In particular, the equations Zx/ + Ljb^b — and 
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Lyi + LibYb = have to hold. Thus, we can express the interior coordmates as 



yi = -L LibYb- 

For non-zero weights uj the matrix Lu is irreducible (that is, the underlying graph is connected, 
see Lemma 3.2) and diagonally dominant. As a consequence L// is invertible and ^ has a unique 



solution ^5, page 363]. 

The barycentric embedding assures that the interior vertices are in equilibrium. However, to 
make G{p) liftable we have to guarantee the equilibrium also for the vertices on /q. We define the 
vectors F :— Fi, . . . ,Ffc as the non-resolving "forces", which arise at the boundary vertices and 
cannot be canceled by the interior stresses: 

Our goal is to define the yet unassigned stresses on the boundary edges such that they cancel the 
forces in F. However, this is not always possible, depending on the shape of the outer face. In 
order to pick a good embedding of /o, we have to know how changing the coordinates of the outer 
face changes the forces in F. The following lemma helps to express this dependence. 

Lemma 3.1 (Substitution Lemma). There are weights Uij = ujji, for i,j G B, independent of the 
location of the boundary vertices, such that 

The weights Qij are the off-diagonal entries of LbiL^^Ljb ^ Lbb- If ^ integral, each Q is a 
multiple of 1/ det L. 

Proof. Let F^, denote the vector (Ff , . . . , Fj^)"^ , where Ff is the x-component of the vector F^. 
We rephrase ([t]) as F^; = Lbb^b + Lbi^i- With help of ^ we eliminate x/ and obtain 

Fa; = Lbb^b ^ LbiL ^Ljb^B ~'- Lxb- 

(The matrix L = Lbb — LbiL^^Ljb is the Schur complement of L in L.) For the y-coordinates, 
we obtain a similar formula with the same matrix L. We define diij as the off-diagonal entries —lij 
of L. Since Lbi = {Ljb)'^ , the matrix L is symmetric and therefore = uiji holds. 

To show that the expression = Lxb has the form stated in the lemma we have to check 
that all row sums in L equal 0. Let 1 denote the vector where all entries are 1, equivalently 
denotes the vector that contains only zeros as entries. Since each of the last n — k rows of L 
sums up to we have LI + Ljb^ — 0; and hence —L^^Lib^ = 1. Plugging this expression into 
LI = Lbb^ ~ LbiL~^Lib^ gives us L\. — Lbb^ + Lbi^, which equals 0. The matrix L can be 
written as a rational expression whose denominator is the determinant of L, and thus the weights 
Co are multiples of 1/ det L. □ 

In linear algebra terms, the lemma can be rephrased as saying that the Schur complement of a 
submatrix of a weighted Laplacian, if it exists, has again the form of a weighted Laplacian. 

The proof assumes that L is invertible. This is the case whenever the graph G has no connected 
component that is a subset of /. If such components exist, they can simply be omitted, since they 
are completely disconnected from B and hence have no effect on the forces in F . Hence, the lemma 
holds for arbitrary graphs, without any connectivity assumptions. 
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Example. For the example of Figure [T] we obtain the following substitution stresses. 

/-(I-12 - ciJia wi2 Wis \ /-96/95 3/5 39/95 \ 

-L = U)i2 -W12-W23 <^23 \^\ '^/^ 

V cZ'ia W23 -Qi3-i:j23j V 39/95 3/5 -96/95/ 

We emphasize that the uj values are independent of the location of x^ and yb- they only depend 
on the combinatorial structure of G. In other words, the stresses (D.y contain all the necessary 
information about the combinatorial structure of G. Thus, we have a compact description (of size 
(2)) of the structure of G that is responsible for the forces in F. We call the stresses uj substitution 
stresses to emphasize that they are used as a substitution for the combinatorial structure of G. 

For the later analysis of the grid size it is necessary to bound the size of the substitution 
stresses. We first state a technical lemma. 

Lemma 3.2. Removing all vertices and edges of a face fo from a 3-connected planar graph G 
leaves a connected graph. 

Proof. After realizing G as a polytope, the claim becomes a special case of the well-known state- 
ment that a graph of a polytope in any dimension remains connected if the vertices of some face 
are removed. This statement can be proved by defining a linear objective function that realizes the 
minimum entirely on the removed vertices. Every remaining vertex is connected to the maximum 
vertex by a monotone increasing path. The objective function can be perturbed such that there is 
a unique maximum. □ 

Lemma 3.3. 1. Let uj be a stress that is positive on every interior edge. Any induced substitu- 
tion stress uJij is positive. 

2. Let Lu be the stress that is 1 on every interior edge. Any induced substitution stress Qij is 
smaller than n ~ k. 

Proof. The substitution stresses are independent of the location of fo- Therefore, we can choose 
the positions for the boundary vertices freely. We place vertex Vi at position (0, 0)-^ and all other 
boundary vertices at (1,0)-^. All vertices he on the segment between (0,0)-^ and (1,0)-^, which is 
the convex hull of the boundary vertices. 

We now show that all interior vertices lie in the interior of this segment. If an interior vertex 
lies at (1,0)"^ then all its neighbors have to lie at (1,0)"^ as well. Otherwise the vertex cannot be in 



equilibrium. But since due to Lemma 3.2 all interior vertices are connected by interior edges this 
would imply that all interior vertices must lie at (1,0)"^. This is a contradiction, since Vi is also 
the neighbor of an interior vertex. By the same arguments one can show that no interior vertex 
can lie at (0, 0)"^. Therefore, all interior vertices have a positive ^-coordinate strictly smaller than 
1. 

In our special embedding the force Fj {j ^ i) can be expressed as Fj — (a)y,0)^. By ^ we 
have UJij = J^k^i '^jk{xj — Xk). Due to the results of the previous paragraph, this sum consists of 
at most |/| summands, which are positive, and in the case a; = 1 smaller than 1. Both statements 
of the lemma follow. □ 

We are now ready to introduce the embedding algorithm. As a first step we construct a 2d 
embedding in equilibrium with respect to a stress w with w = 1 on the interior edges. In order to get 
equilibrium on the boundary vertices as well, we have to choose their locations and the stresses on 
the boundary edges appropriately. This leads to a non-linear system in the 2k unknowns x^, and 
Yb and the k unknown boundary stresses 0^12, a;23, . . . ,a;fci. Let Lq be the Laplacian of the graph 
that consists of the outer face /g only, with unknown stresses a;i2, W237 • ■ • jWfei for the boundary 
edges. The 2k equations of the system are given by 

LqXb + LxB = 0, Loys + iys = 0. (8) 
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Since these equations are dependent, the system is under-constrained. To solve it, we fix as many 
boundary coordinates as necessary to obtain a unique solution. We also have to ensure that the 
solution defines a convex face. We continue with a case distinction on k. 



Case 1: G contains a triangular face 

The triangular case is easy: we can position the boundary vertices at any convenient position (see 
for example [U). We choose: 




Lemma 3.4. // G contains a triangle and we place the boundary vertices as stated in ([9]) then the 
boundary forces can be resolved. 

Proof. We embed G as barycentric embedding and calculate the substitution stresses. After setting 
L0i2 = — a)i2,W23 — —0)23, Wis = — all points are in equilibrium. □ 



Case 2: G contains a quadrilateral but no triangular face 

If /o is a quadrilateral we have to fix some coordinates of the boundary vertices such that it is 
possible to cancel the forces in F. We used computer algebra software to experiment with various 
possibilities to constrain the coordinates and solve the non- linear system ([S]). A unique solution 
can be obtained by setting 

Pi = (0) , P2 = (J) , P3 = {1)1 , P4 = (J) , (10) 



with 



2/3 = TP —■ (11) 

The solution of the equation system ([8| provides also the stresses on the boundary edges. These 
stresses are not necessary for our further computations; we mention them here for completeness 
only. 

LO12 = — 2wi3 — a)i2, 

CJ23 = '^24 ^ 2ci)i3 — 0)23, 

^24 - ^12) 

cj34 = -^ — W34, y-^) 

'^24<^13 
Wl4 — TT^ W14. 

a;24 — ^1^13 

We assume that 0)13 > tli24. (Otherwise we cyclically relabel the vertices on /q.) Since w = 1 on 
the interior edges the substitution stresses are positive by Lemma |3.3[ Under this assumption we 
can deduce that < 1/3 < 1- Hence, /o forms a convex face. 

Note that the substitution stresses Uij between adjacent vertices (on the boundary) are irrel- 
evant. The forces resulting by the boundary stresses uJij can be directly canceled by the corre- 
sponding stresses uiij. This can also be observed by looking at the solution of the corresponding 
equation system: Boundary stresses do not appear in the solution for j/3. (Furthermore the sum 
Cbij + UJij for boundary edges (i,j) does not depend on any other boundary stress either.) 



Lemma 3.5. If G contains a quadrilateral and we place the boundary vertices as stated in ( 10 ) an 



(111, then /o forms a convex quadrilateral and the boundary stresses (12 I cancel the forces F. □ 
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Case 3: G contains no triangular and no quadrilateral face 

The case if the smaUest face of G is a pentagon is more complicated. We have (2) = 10 substitution 
stresses w^j, but the adjacent ones do not count (by the same reasons given in the previous case). 
So we are left with five "diagonal" substitution stresses W13, (2;i4, (I;24, and cl^as. 

Like in the previous cases we determine a unique solution of the equation system by fixing 
some of the coordinates of the outer face. However, we have to make more effort to guarantee the 
convexity of /q. We first observe: 

Lemma 3.6. We can relabel the boundary points for any stress (wy)i<ij<5 such that 

'^35 > '^24 and UJ25 > ^Id- 
Proof. Without loss of generality we assume that the largest stress on an interior edge is lu^c^. If 
W25 > 1^13 we are done. Otherwise we relabel the vertices by exchanging p3 ps and Pi P2- □ 

For the rest of this section we label the vertices such that Lemma [3l6] holds . The way we embed 
/o depends on the substitution stresses oJij . 
Case 3A: 

We assume that 



'^35'^14 + 1^14W25 + W25W24 + "^13^35 > 



(13) 



In this case we assign 
Pi 



,P2 



,P3 



,P4 



,P5 



Figure [2]Ja) illustrates the location of the points. Together with the equations of ([8| we obtain as 



P4 



P3 



Pi 



P2 



P5 



Pa 



P3 



P2 



(a) Case 3A 

Figure 2: Placement of the boundary vertices 



Pi 

(b) Case 3B 



solution for ps: 



X5 



(tills - ^25 - t^24)(t^35 + ^^13 - ^^24) 
W35'^14 + ^li^25 + 'I'25'^24 + Wl3'^35 ~ '^'35^25 ' 
W35 + W13 — OJ24 



2/5 



W35 + W25 



The boundary stresses cjij are complicated expressions that are not necessary for further compu- 
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tations. We list them here for completeness only. 

'^13('^i5 + ^^24^35 + 20)24^25 " 'il3'^25) +tl'l4(w|5 + W25'^35 + ^24^25 + W35W24) 

W12 ~ z z z z ; z z z z z W12, 

"^35^25 ~ '^14<^25 — W25W24 — W13CLI35 — a;35a;i4 

Wl4(wi5 + ^^35^13 + W25'^35 + '^13W25)+W24(wi5 + ^130)25 + 2a;i3W35 - W35W24) 
W34 — z ; : : z z z z z z OJ34, 

W35a;25 — Wi4a;25 — W25W24 — — UJ35LV14 



Wi3a;25 + W25a;35 + W24W25 

W23 — z : 1^23 7 

—'^25 ~ 1^35 

W24t<j25 + W25W14 + (I'i4ii'35 + '^24^^35 

W45 = z ; z W45 , 

W13 — a;24 — W25 



Wi3a;25 + <^35Wl3 + Wl4'^25 + '^14'^35 

Wi5 — z z z Wl5- 

^24 — W35 — W13 

We have to check that /o forms a convex polygon. Clearly, > 0, since the w^ 's are greater 
than zero and 0)35 > a;24. Moreover < 1, because UJ25 > ^13. The numerator of X5 is negative 



and due to (13 1 the denominator of X5 is positive. Therefore, 3:5 < and /o forms a convex 
polygon. 

Case 3B: 



We assume the opposite of (131, namely 

t^35t^i4 + W14W25 + W25W24 + '^i3(I;35 < a)35'^25- (14) 

The coordinates for the boundary vertices are chosen as 

0\ /l\ /l\ fO\ (-1 



See Figure [2jb) for an illustration. This leads to the solution 



2/2 = -2 



tl)24(l'l3 + £^24^^35 + '^25'^13 + 2tl)25'^35 " '^{^ - 2wi3<^35 - '^Z'o^ 
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W24'^35 + (i'25'^13 + 2c2;25'^35 

i i)'^ aI n -A- { I In A f I I'D r:: - 

2/3= 2 



'^24^^13 + t<j24W35 + (I'25Wl3 + 2(I)25t2'35 — UJ24 ~ 20)24^25 ^ 'il4'i25 



W24W35 + (^25^13 + 2a;25'^35 

The boundary stresses LUij are once more not necessary for further computations and listed for 
completeness only. 

W12 = — W24 — 2(1)25 — W12, 

-W25('^?3 + 2wi3'^35 + 2(Ii24tI'35) " l^S^if^h + £^25^14) 
2^35 (£l;24 + ^25 -^^13 - I^^m) + 2^25 (^^13 +£^35-^^24 " ^'^14) - (^13 - ^24)^ 

1^34 ~ —'^14 ~ 2a)i5 — W34, 

^45 = 'i'24 ~ 2W35 — a)i3 — UJ^^, 

^15 = 'i'13 ~ 2W25 ~ ^1)24 — (^15- 

The outer face is convex if —2 < 7/2 < 2/3 < 2. The inequalities — 2 < 7/2 and 2/3 < 2 are 
equivalent to 

-0)13 - W35W14 + 0)13(^24 - 20)35) < and 
-0)24 ~ '^14^25 + W24('^i3 - 20)25) < 0. 



10 



Both inequahties hold, because we add only negative summands on the left side. It remains to 
check if y2 — Us < 0. First we get rid of the denominator and bring all negative summands on the 
right side. This leads to the equivalent inequality 

'^IS + '^24 + 2W13'^35 + 2(I'24'i25 + ^^25^14 + WasW^ < 

2^24^35 + 2a)25'^l3 + 4,1^25^35 + 20)24^13 • (15) 



We observe that uj^^ < W25W13 and uj2i ^ W24W35. Because of the assumption (14 1 we have 



4£l'35a'25 > 20)130)35 + 2ui24U)25 + W25W14 +^350)14. Therefore, the right side of (15) is greater than 
its left side, which shows that 7/2 < 2/3 and fo forms a convex pentagon. This completes the case 
distinction and we conclude with the following lemma. 

Lemma 3.7. // we place the boundary vertices as discussed above, then the outer face will be 
embedded as a convex pentagon and the computed boundary stresses cancel the forces in F. □ 

We have defined four different ways to embed G. The selected embedding depends on the 
combinatorial structure G. If G contains a triangular face we say it is of type 3. If it contains 
a quadrilateral but no triangular face G is of type 4- Otherwise the embedding depends on the 



substitution stresses induced by the combinatorial structure of G. If ( 13 1 holds (Case 3A) G is of 
type 5A, otherwise (Case 3B) we say G is of type 5B. 

3.2 Lifting the Plane Embedding and Scaling to Integrality 

We continue with lifting the plane embedding of G to M-^. With help of the observations made 
in Section [2] the incremental computation of the 3d embedding is straightforward. It suffices to 
compute for every face fi the corresponding plane Hi. 

Since the embedding of /o is convex, the boundary stresses must necessarily be negative, since 
otherwise the boundary vertices could not be in equilibrium with all interior stresses being positive. 
Thus we do not need to explicitly check the sign of the boundary stresses. The sign pattern of the 
stress implies that the lifting of the plane embedding gives a convex polytope (see Proposition [T]) . 

As described in Section [2] we begin the lifting by fixing the plane Hi for some interior face /i 
as the x-y-plane. We set ai = (0, 0)^, di = 0, and compute the remaining planes face by face using 
equations ^ and ([5]) . It is not necessary to compute the parameters of TJq since we can determine 
the heights of pi, . . . , pfc by some plane Hi of an interior face. Hence, the lifting can be computed 
using only stresses on interior edges. This simplifies the later analysis because all interior stresses 
are 1, whereas the boundary stresses are complicated expressions. 

It can be observed that the computed lifting has rational coordinates. This is true because the 
barycentric embedding gives rational coordinates and the lifting process is based on multiplication 
and addition of the 2d coordinates. Hence, the z-coordinates are also rational. We analyze the 
common denominator of the coordinates to obtain scaling factors for the integral embedding. We 
use different scaling factors Sx for the x-coordinates and Sy for the y-coordinates. 



As a consequence of Lemma 2.1 it is sufficient to scale to integer x and y-coordinates. Further- 
more we observe: 

Lemma 3.8. // the boundary points are integral, the barycentric embedding yields coordinates that 
are multiples of 1/ det L. 

Proof. The interior plane coordinates are a result of equation ([6|. By Cramer's rule every coordi- 
nate can be expressed as 

x^ = det / det L, 

where detL^*^ is obtained from L by replacing the i-th column of L by Ljb^b- Since detZ^*' is 
integral, detZ is the denominator of Xi. The same holds for yi. □ 
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Our first goal is to scale the plane embedding such that the boundary vertices get integer 
coordinates. Let be the integral scaling factor that gives integer boundary x-coordinates and 
Sy be the integral scaling factor that gives integer boundary y-coordinates. Due to Lemma 



3.8 



the scaling factors Sx ■— det L and Sy ■= Sy det L make all vertices integral. Since we choose 
integral scaling factors S^ and Sy no integer coordinate is scaled to a non-integer. 

Let us now compute the factors that are necessary to scale to integer boundary coordinates. 
Clearly the scaling factors depend on the type of G. If G is of type 3 then we need not scale, since 
all boundary coordinates are either or 1. If G is of type 4 we have to scale the ^-coordinates only 
(see |ll|). We multiply with Sy := (2a;i3 — ^^24) det L, which gives S'^ya = W24 det L, which due 
to the substitution Lemma is an integer. 

If G is of type 5A we have to scale such that X5 and 2/5 become integral. We pick 

= (^35^14 + £^14^25 + ^25^^24 + ^laWas - a)35W25)(det Z)^, 

Sy = (0)35 +cl'25) det L. 

It can be easily checked that these factors as well as S^x^ and S'^ys are integral. 

When G is of type 5B, the only non-integer boundary coordinates are j/2 and 1/3, we need to 
scale in y-direction only. We choose 

Sy = (W24W35 + W25W13 + 2(Ii25(I'35)(detL)^. 

Again, due to the Substitution Lemma, Sy , 6*^2/2, and 5^2/3 are all integral. 

For every type of G there is a pair of scaling factors Sx,Sy, such that the scaled boundary 
points are integral. Table [l] summarizes the discussion and lists the final scaling factors depending 
on the type of G. 



type of G 


scaling factors 


3 


Sx = Sy = det L 




Sx — det L 


4 


Sy = (2^13 -ti24){detZ)2 


5A 


Sx = (^^35^^14 + W14W25 + ^^25^24 + '^13'^35 " a;35a;25)(det Z)^ 

Sy = (0)35 +W25)(detL)2 


bB 


Sx = det L 


Sy = (W24'^35 + W25W13 + 2tJ25W35)(det L)^ 



Table 1: The scaling factors Sx and Sy for the different types of G. 



3.3 Analysis of the Grid Size 

To bound the size of the coordinates of the integer embedding it is crucial to obtain a good bound 
for detL. Recall that we assume unit stresses w = 1 on the interior edges, throughout. There 
exists a connection between the number of spanning trees in G and det L. Let us first define: 

Definition 2. Let B be a subset of vertices of G. A subgraph of G is called spanning S-forest if 

• it consists of \B\ vertex disjoint trees covering all vertices of G, 

• each tree contains a unique vertex from B. 
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In the following we use the set of boundary vertices for B. Let Tb{G) denote the number of 
spanning _B-forests of G and T{G) the number of spanning trees of G. A generalization of the 
Matrix- Tree Theorem jT6j (see also [21]) states that the number of spanning _B-forests of G is det L. 
In our case, we can directly bound Fb{G) by T{G). 

Lemma 3.9. Let G he a planar graph with a distinguished face and let B he the set of vertices of 
this face. The numher of spanning B-forests of G is hounded from above hy 

J'siG) < T[G). 

Proof. Every spanning i3-forest can be turned into a spanning tree by adding all boundary edges 
except (1,2). No two distinct spanning _B-forests are associated with the same spanning tree. 
Therefore, the number of spanning trees exceeds the number of spanning _B-forests. Since there is 
a spanning tree that contains the edge (1, 2) the inequality is strict. □ 

It is easy to give an exponential upper bound for T{G): 

Proposition 2 (Ribo Mor [21] ). 1. The numher of spanning trees in a graph is hounded hy the 
product of all vertex degrees: 

nG)<\{deg{v.). 

i 

2. For a planar graph, we have T{G) < Y[i deg(ui) < 6". 

Proof. 1. Consider all directed graphs that are obtained by choosing an outgoing edge in G out of 
every vertex except u„. The number of these directed graphs is given by rir=i^ deg(ui). By ignoring 
the edge orientations, one obtains all spanning trees (and many graphs that are not spanning trees). 
Alternatively, the bound can be proved by applying a variant of Hadamard's inequality for positive 
semidefinitc matrices [ST, page 176] to the (positive scmidcfinitc) matrix V that is obtained by 
removing from the Laplacian L the row and column corresponding to the vertex v^'. 

n—l n—1 

T{G) = detZ < n = Jl deg(w,) 

i=l i=l 

2. This follows from the arithmetic-geometric-mean inequality and the fact that ^ • deg(wi) < 
6n, which is a consequence of Euler's formula. □ 

Sharper bounds for T(G) have been given by Ribo Mor et al. [33], see also [3T1 [51]. These 
bounds take into account whether G contains triangular or quadrilateral faces: 



if G is of type 3 
if G is of type 4 
if G is of type 5A/5B 



Tb{G) < T{G) < 5.3", 
Tb{G) < T{G) < 3.529988'^ 
Tb{G) < T{G) < 2.847263'^ 



Since we know upper bounds for the a) values (by Lemma 3.31 and detX (by the previous 
discussion) we can bound the size of the integer coordinates of the embedding of G. We start with 
bounding the x and y-coordinates. Let Ax denote an upper bound for the difference between the 
largest and the smallest cc-coordinate. Ay is defined in the same way for the y-coordinates. 

Again we have to discuss the 4 cases separately. If G is of type 3 then clearly Ax = Ay = 
det L. If G is of type 4 the largest x-coordinate is 2Sx and the smallest zero. Thus we have 
Ax = 2detZ. The largest y-coordinate is obtained at y4 = 1 (remember < 1), therefore 
Ay ~ Sy ~ (2(Iii3 — a)24)(det L)^. Let us now assume G is of type 5A. The value of Ax is given by 
X2 — X5. Evaluating this expression leads to 

Ax = (W25(wi3 +W14) +UJ35{uJi4 -I- (1)25) " (^^13 - W24)^)(det Z)^. 
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type of G Ax 



3 detL 

4 2 det L 

5A (aj25('^i3 + i^id + '^35(0)14 + W25) - (wi3 - a;24)^)(det L)^ 

5B 2 det L 



Table 2: The values Ax depending on the type of G. 



type of G Ay 

3 detZ 

4 (2W13 - W24)(detZ)2 

5A (^35 +W25)(detZ)2 

5B 4(W24W35 + W25'^13 + 2aj25W35)(det L)^ 

Table 3: The values Ay depending on the type of G. 

Since the smallest ^/-coordinate is zero we have Ay — y^, which equals + 0)25 ) (det L)^. It 
remains to discuss the case when G is of type 5B. Before the scaling the coordinates fulfill 
— 1 < x < 1 and —2 < y < 2. Combining these inequalities with the scaling factors yields 
Ax — 2 det L and Ay = 4,{uj24^35 + '^25^^13 + 2(Ii25W35)(det L)^ . We sum up the results for Ax and 
Ay in Table [2] and Table [3] With help of Lemma [331 we can eliminate the oj values that appear in 
the bounds of Ax and Ay. The resulting upper bounds, which we use in the further analysis, are 
listed in Table H 

We finish the analysis of the necessary grid size by calculating the size of the z-coordinates. 

Lemma 3.10. Let G(p) be an integral 2d embedding of a graph with n vertices with equilibrium 
stress uj and let the stress on all interior edges be 1. The difference between two x- coordinates is 
at most Ax and the difference between two y -coordinates is at most Ay. Then we have an integral 
lifting with 

< Zi < 2nAxAy 

for all z-coordinates Zi. 

Proof. Due to Lemma |2.1| we know that there exists an integral lifting for the setting described 
in the lemma. We place an interior face fi in the xy-plane and compute the lifting by using the 
stresses on the interior edges. Notice that all z-coordinates are non-positive in this lifting. Thus it 
suffices to compute the smallest z-coordinate. The claimed lifting is then obtained by translating 
the polytope such that the smallest z-coordinate becomes 0. 

We choose as face /i a face that shares an edge with the outer face /q. Furthermore we assume 
that the boundary point farthest away from the line that contains fi D /o is located in the origin 



type of G upper bound for Ax upper bound for Ay 

3 det L det L 

4 2detL 2n(detZ)2 
5 A 4n2(detL)3 2n(detZ)2 
5B 2detL l6n'^{detLf 

Table 4: Upper bounds for Aa; and Ay depending on the type of G. 
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(let this point be pi). This is no restriction since a translation of the embedding does not interfere 
with the lifting. The lifted polytope lies below the a;y-plane Hi and above Hq. We notice that the 
smallest z-coordinatc of Hq (and hence the smallest z-coordinatc of the embedding) is realized at 
Pi 

Let fk be an interior face that contains pi. The z-coordinate of pi is given as 

zi = (ak,pi) +dk = dk- 

The variable dt can be computed with help of equation Let C be a set of interior edges 
that are crossed by "walking" from /i to fk- Due to Euler's formula G has at most 2n — 4 faces. 
No face is entered twice and thus every face contributes at most one edge (the edge where the 
"walk" leaves the face) to the set C. This implies that C includes at most 2n — 3 edges. We ignore 
the orientation of the edges at this place since it does not matter for bounding dk- We deduce 

~dk< |(p,,PjL)| < 2nmax{|(p„p^-L)|: 1 < i,j < 7i}. 

For two points pi , pj we have 

Thus, (pi, pj-) equals two times the negative area of the triangle spanned by p^, pj, and the origin 
(which coincides with pi). This triangle is contained inside the embedded outer face /o and also 
inside a rectangle with edge lengths Ax and Ay. A rectangle has at least twice the area of an 
inscribed triangle. To see this, observe that an inscribed triangle with the largest area must have 
one of the rectangle edges as base and the other as height. Thus |(Pi,Pj")| < Aa;Aj/ and the 
smallest z-coordinate is larger than — 2nAxAi/. □ 

By applying Lemma |3.10[ we compute the bounds for the z-coordinates, using the values of 
Aa; and Ay listed in Table |4] We conclude with the main theorems. 

Theorem 3.1. Every 3-polytope with n vertices whose graph contains at least a triangle can be 
realized on an integer grid with 

< x,,y, < 5.3", 
0< z, <2n-28.4". 

Theorem 3.2. Every 3-polytope with n vertices whose graph contains at least one quadrilateral 
face can be realized on an integer grid with 

< X, < 2 • 3.530", 
0<y,<2n- 12.461", 
< z, < • 43.987". 

For the most general theorem we have to combine the analysis for the cases 5 A and 5B. We 
rotate the embedding if G is of type 5B by exchanging its x and y-coordinates to obtain a better 
bound. The largest z-coordinate is given by max{167i'^ • 187.128", 64^3 • 65.722"} = 1671-* • 187.128", 
since 7i > 4 if G contains a pentagon. Thus the largest bound on the z-coordinate arises from 
case 5^. 

Theorem 3.3. Every 3-polytope with n vertices can be realized on an integer grid with 

0<Xi < 16n^ ■ 23.083", 

< yi < 2n • 8.107", 

< Zi < 16n'^ ■ 187.128". 
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We can improve the constant factor of the z-coordinate by a more careful analysis. This can be 
achieved by placing the face fo in the a;y-plane and then compute the lifting using the interior edges 
but also one boundary edge. As mentioned before the structure of the stresses on the boundary 
edges is more complicated. Since the improvement would only be a constant factor we decided to 
present the easier analysis with help of Lemma [3.1 0| The more complicated analysis can be found 
in [12]. 

We see that exponentially large coordinates suffice to embed G as 3-polytope. The exponential 
growth of the size of the coordinates is determined by (detZ)''^. 

Corollary 3.1. Every 3-polytope with n vertices can be realized with integer coordinates of size 

Let us add some remarks on the running time of the embedding algorithm. If we know the 
substitution stresses the computation of the location of the outer face can be done in constant time. 
The same is true for the scaling factors. Once we computed the plane embedding, the lifting can 
be computed face by face, which needs in total 0(n) steps. The computation of the substitution 
stresses and of the interior vertices can be done by solving a linear system. Since its underlying 
structure is planar, we can use nested dissections based on the planar separator theorem to solve 
it [TTIHH]- This implies that a solution can be computed in 0{M{^/n)) time, where M(n) is the 
upper bound for multiplying two nxn matrices. The current record for M(n) is 0(n^-^^^) which is 
due Coppersmith and Winograd [SI. Thus the overall running time is given by 0{v}'^^^) arithmetic 
operations. 

4 An Example: the Dodecahedron 

The regular dodecahedron is one of the five Platonic solids. It has 20 vertices, 30 edges and 12 
faces, which are regular pentagons. Figure [3| shows the graph and a 3-dimensional realization of 
it. It is the smallest polytope without triangles and quadrilateral faces, and thus we have to apply 
the more involved cases. Since the dodecahedron is symmetric it makes no difference which face 




(a) (b) (c) 

Figure 3: (a) The graph of the dodecahedron and (b,c) two realizations as 3-polytope. 

we choose as the outer face. We start the computation with calculating the uj values. Remember, 
these values are the off-diagonal entries of the matrix —{Lbb — LbiL~^Lib)- We obtain for all the 
stresses W13, <Di4, 0)24, ^^25 a-nd 0)35 the value 36/449. The fact that all these stresses have the same 
value is again due to the symmetry of the dodecahedron. Since the outer face is a 5-gon, G is of 
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type 5A or 5B. Evaluating (13) shows that the graph is of type 5A. With help of the substitution 
stresses we compute the coordinates of the boundary vertices. We obtain 




Figure 4: The barycentric (plane) embedding of the dodecahedron. 



This yields the scaling factors 

^ 1 264 158 727 403 904, Sy = 26 069 428 512. 

We continue with the lifting of the plane embedding to M'^. The faces are lifted incrementally 
as described in Section [2j The numeric data of the lifting are listed in [29] . Figure [5] shows the 
computed embedding. We have scaled down the z-coordinates to obtain an illustrative picture. 
The highest absolute coordinate is 



I Z3 1 = 1 1 083 163 098 782 678 334 820 352 « 2 



83.19 



which is smaller than the bound 2^^^ from Corollary 3.1 

The computed embedding allows a smaller integer realization. Due to the fact that the greatest 
common divisor of the x-coordinates is 938 499 426 432 = 449^ x 10365 and the greatest common 
divisor of the y-coordinates is 29 030 544 — 449^ x 144, scaling down by these factors yields a smaller 
integer embedding. We obtain an integral plane embedding on the grid [—27, 1347] x [0, 898]. The 
corresponding z-coordinates range between and 406 497. This reduction is due to the fact that all 
substitution stresses w are equal. Thus one might have replaced them by a; = 1 in the subsequent 
calculations. 

A much smaller grid embedding of the dodecahedron was constructed by hand by Francisco 
Santos. It is centrally symmetric and fits inside a 6 x 4 x 8 box, see Figure [3][b). It is hard to 
believe that a smaller realization would be possible. Another, more symmetric, realization of the 
dodecahedron is the pyritohedron (one of the possible crystal shapes of the mineral pyrite), as 
pointed out to us by Gabor Gevay. It fits in a 12 x 12 x 12 box, see Figure [3]jc). It has 8 vertices of 
the form (±4, ±4, ±4), plus 12 vertices, which are the 4 vertices of the form (0,±3, ±6) and their 
cyclic rotations of the coordinates. The normals of the 12 faces are the vectors (0,±2, ±1) and 
their cyclic rotations. 
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Figure 5: Two views of the dodecahedron embedded with our algorithm, with scaled 2;-axis. The 
right picture includes also the equilibrium-stressed plane embedding. 

Acknowledgements. We thank a referee for a very thorough reading of the manuscript. 

References 

[1] D. M. Ackcta and J. D. Zum'c. On the maximal number of edges of convex digital polygons 
included into an 'm x m-grid. J. Comb. Theory Ser. A, 69(2):358-368, 1995. 

[2] G. E. Andrews. A lower bound for the volume of strictly convex bodies with many boundary 
lattice points. Trans. Amer. Math. Sac, 99:272-277, 1961. 

[3] I. Barany and G. Rote. Strictly convex drawings of planar graphs. Documenta Math., 11:369- 
391, 2006. 

[4] N. Bonichon, S. Felsner, and M. Mosbah. Convex drawings of 3-connected planar graphs. 
Algonthmica, 47:399-420, 2007. 

[5] K. Buchin and A. Schulz. On the mimbcr of spanning trees a planar graph can have. In 
Proceedings of the 18th Annual European Symposium on Algorithms, (ESA 2010), volume to 
appear, 2010. 

[6] M. Chrobak, M. T. Goodrich, and R. Tamassia. Convex drawings of graphs in two and three 
dimensions (preliminary version). In 12th Annual Symposium on Computational Geometry, 
pages 319-328, 1996. 

[7] R. Connelly, E. D. Demaine, and G. Rote. Straightening polygonal arcs and convexifying 
polygonal cycles. Discrete and Computational Geometry, 30:205-239, 2003. 

[8] D. Coppersmith and S. Winograd. Matrix multiplication via arithmetic progressions. J. Symb. 
Comput., 9(3):251-280, 1990. 

[9] H. Crapo and W. Whiteley. Plane self stresses and projected polyhedra I: The basic pattern. 
Structural Topology, 20:55-78, 1993. 

[10] G. Das and M. T. Goodrich. On the complexity of optimization problems for 3-dimensional 
convex polyhedra and decision trees. Comput. Geom. Theory AppL, 8(3): 123-137, 1997. 



18 



[11] H. dc Fraysscix, J. Pach, and R. Pollack. How to draw a planar graph on a grid. Combinatorica, 
10(1):41 51, 1990. 

[12] P. Eades and P. Garvan. Drawing stressed planar graphs in three dimensions. In F.-J. 
Brandenburg, editor, Graph Drawing, volume 1027 of Lecture Notes in Computer Science, 
pages 212-223. Springer, 1995. 

[13] S. J. Gortlcr, C. Gotsman, and D. Thurston. Discrete one-forms on meshes and applications 
to 3d mesh parameterization. Computer Aided Geometric Design, 23(2):83-112, 2006. 

[14] J. E. Hopcroft and P. J. Kahn. A paradigm for robust geometric algorithms. Algorithmica, 
7(4):339-380, 1992. 

[15] R. A. Horn and C. R. Johnson. Matrix Analysis. Cambridge University Press, 1990. 

[16] M. Lewin. A generalization of the matrix-tree theorem. Mathematische Zeitschrift, 181(1):55- 
70, 1982. 

[17] R. J. Lipton, D. Rose, and R. Tarjan. Generalized nested dissection. SIAM J. Numer. Anal., 
16(2):346-358, 1979. 

[18] R. J. Lipton and R. E. Tarjan. Applications of a planar separator theorem. SIAM J. Comput., 
9(3):615-627, 1980. 

[19] J. C. Maxwell. On reciprocal figures and diagrams of forces. Phil. Mag. Ser., 27:250-261, 
1864. 

[20] S. Onn and B. Sturmfels. A quantitative Steinitz' theorem. In Beitrdge zur Algebra und 
Geometric, volume 35, pages 125-129, 1994. 

[21] A. Ribo Mor. Realization and Counting Problems for Planar Structures: Trees and Linkages, 
Polytopes and Polyominoes. PhD thesis, Preie Universitat Berlin, 2006. 

[22] A. Ribo Mor, G. Rote, and A. Schulz. Embedding 3-polytopes on a small grid. In SCG'07: 
Proc. 23rd Annual Symposium on Computational Geometry, pages 112-118, New York, NY, 
USA, 2007. ACM. 

[23] A. Ribo Mor, G. Rote, and X. Yong. Upper bounds for the number of spanning trees of a 
planar graph, in preparation, 2009. 

[24] J. Richter-Gebert. Realization Spaces of Polytopes, volume 1643 of Lecture Notes in Mathe- 
matics. Springer, 1996. 

[25] J. Richter-Gebert and G. M. Ziegler. Realization spaces of 4-polytopes are universal. Bull. 
Amer. Math. Soc, 32:403, 1995. 

[26] G. Rote. The number of spanning trees in a planar graph. In Oberwolfach Reports, volume 2, 
pages 969 973. European Mathematical Society - Publishing House, 2005. 

[27] W. Schnyder. Embedding planar graphs on the grid. In Proc. 1st ACM-SIAM Sympos. 

Discrete Algorithms, pages 138 148, 1990. 

[28] O. Schramm. Existence and uniqueness of packings with specified combinatorics. Israel J. 
Math., 73:321-341, 1991. 

[29] A. Schulz. Lifting Planar Graphs to Realize Integral 3-Polytopes and Topics in Pseudo- 
Triangulations. PhD thesis, Preie Universitat Berlin, 2008. 



19 



[30] A. Schulz. Drawing 3-polytopcs with good vertex resolution. In D. Eppstein and E. R. 
Gansner, editors, Graph Drawing, volume 5849 of Lecture Notes in Computer Science, pages 
33-44. Springer, 2009. 

[31] E. Steinitz. Polyeder und Raumeinteilungen. In Encyclopddie der mathematischen Wis- 
senschaften, volume 3-1-2 (Geometrie), chapter 12, pages 1-139. B. G. Teubner, Leipzig, 
1916. 

[32] T. Thiele. Extremalprobleme fiir Punktmengen. Master's thesis, Freie Universitat Berlin, 
1991. 

[33] W. T. Tutte. Convex representations of graphs. Proceedings London Mathematical Society, 
10(38) :304-320, 1960. 

[34] W. T. Tutte. How to draw a graph. Proceedings London Mathematical Society, 13(52) :743-768, 
1963. 

[35] K. Voss and R. Klcttc. On the maximal number of edges of convex digital polygons included 
into a square. Pocitace a umeld inteligencia, l(6):549-558, 1982. (in Russian). 

[36] W. Whiteley. Motion and stresses of projected polyhedra. Structural Topology, 7:13-38, 1982. 

[37] F. Zhang. Matrix Theory. Springer, 1999. 

[38] F. Zickfcld. Geometric and Gombinatorial Structures on Graphs. PhD thesis. Technical 
University Berlin, December 2007. 



20 



